Carbon sequestration characteristics of two plantation forest ecosystems with different lithologies of karst

In karst regions, the majority of studies have focused on ecosystem carbon sequestration in the same lithology, but no studies in different lithologies. In this study, actual measurements were used to reveal carbon sequestration characteristics of two plantation forest ecosystems (Bodinieri cinnamon and Cupressus funebris) with different lithologies of karst. The results showed that the tree layer showed the highest vegetation biomass, carbon content, carbon density, and ratio of aboveground biomass to belowground biomass. The carbon density of B. cinnamon plantation and C. funebris plantation was high in dolomite and in limestone respectively. The soil quality and carbon density of bare ground and plantation varied across different lithologies. The carbon density of various ecosystem components was in the order of vegetation>soil>litterfall. The carbon density and net carbon density of plantation varied across different lithologies. In B. cinnamon plantation, the carbon sequestration rate of vegetation and ecosystem was high in dolomite, moderate in limestone, and low in dolomitic sandstone. In Cupressus funebris plantation, the carbon sequestration rate was in the order of limestone>dolomite>dolomitic sandstone. These findings revealed that lithology is an important factor affecting ecosystem carbon pools, and plantation ecosystems have low biomass and low carbon density in karst areas.


Introduction
Global warming is an unavoidable problem. Scientists are actively exploring carbon peaking and carbon neutral strategies. To achieve carbon peaking and carbon neutrality, it is inevitable to identify the carbon sequestration rate of various ecosystems and their uncertainty, stability, and sustainability. However, these characteristics remain unclear for many ecosystems, particularly the karst regions, due mainly to high landscape heterogeneity, evident dualistic structure, small environmental capacity, discontinuous and shallow soil layer, exposed rocks, poor water storage, and weak resistance to disturbance [1,2]. These challenges further complicate the assessment of carbon sinks in the ecosystems of karst regions, which account for approximately 15% of the global land area and provide drinking water to nearly one-fourth of the global population [3,4]. Therefore, the evaluation of carbon sequestration in karst ecosystems is of great theoretical and practical significance to assess carbon sinks in karst forests, establish carbon sink forests, and implement carbon neutral strategies.
Previous studies related to carbon sinks in the forest ecosystems of karst regions mainly focused on ecosystem carbon storage and its distribution in plantation forests of different ages [5][6][7]; carbon sequestration in urban forests [8,9]; ecosystem carbon sequestration in the natural recovery of degraded karst forest vegetation [10]; carbon fluxes at different stages of the secondary succession of karst grasslands [6]; and carbon sinks of ecosystem components, including vegetation [11], soil [12][13][14][15][16][17], and litterfall [18]; however, there have been no studies on carbon sequestration in plantation forest ecosystems under different lithologies in karst regions. In karst regions, lithology is an important factor driving differences in plant habitats [2,19], and the physicochemical properties and fertility of soil produced through weathering also vary under different lithologies, thus producing a greater impact on plant physiological and ecological functions [20,21]. For instance, karst regions in the Guizhou Province of China, which is at the center of the karst regions of southern China, account for 25.8% of the total karst area in the country, with exposed carbonate rocks covering an area of 13 × 104 km 2 . The Guizhou Province is home to the largest and the most developed tropical and subtropical karst in the world, including China [22]. According to the distribution of carbonate rock fabric [23], there are seven types of lithologies. Among these, limestone, dolomite, and dolomitic sandstone account for 42.73%, 17.31%, and 6.59%, respectively, of the total area of carbonate rocks in the Guizhou Province. Cupressus funebris and Bodinieri cinnamon are the dominant plant species growing in soils developed from these lithologies. Therefore, in the present study, we selected C. funebris and B. cinnamon planted on dolomite, limestone, and dolomitic sandstone in the Guizhou Province and obtained field measurements for studying the carbon sequestration of these two plantation ecosystems to reveal the composition of and changes in plantation ecosystem biomass and carbon density, soil carbon pool characteristics, litterfall carbon pool characteristics, and ecosystem carbon sequestration capacity under different lithologies in karst regions. Our findings can serve as the foundation for carbon sink assessment, ecosystem service function evaluation, carbon sink forest establishment, and vegetation restoration in karst regions.

Overview of the study area
The Guizhou Province is located at the eastern end of the Yunnan-Guizhou Plateau in southwestern China (24˚37 0 to 29˚13 0 N, 103˚36 0 to 109˚36 0 E), with terrain descending in steps from the northwest to southeast and an average elevation of 1,100 m. In this region, karst landforms are typically developed and widely distributed, with distinct morphological types and obvious geographical differentiation, constituting a special karst ecosystem. Therefore, the Guizhou Province is an ideal area for studying karst forest ecosystem restoration. The province has monsoon-influenced humid subtropical climate, and owing to the mountainous terrain, the climate is complex and diverse. The western part belongs to the warm temperate zone, the southern and northwestern parts belongs to the southern subtropical zone, and the rest of the region gradually transitions from the central subtropical zone to the northern subtropical zone, with the raising of Miaoling, Wuling Mountains, and Dalou Mountains. The mean annual temperature and precipitation are 10-18˚C and 1,100-1,500 mm, respectively. The relative humidity is over 70%, and the annual sunshine hours are 1,300. The frost-free period is constant weight in an oven at 80˚C. The ratio of dry mass to fresh mass (kg) was separately calculated for each part. Then, the sample dry mass (kg) was calculated, and the biomass (t�hm -2 ) per hectare of each sample in every sample plot was estimated. The biomass of each harvested organ of the standard trees was subjected to regression analysis using, and the following regression equation was obtained [10,[27][28][29][30]: where W is the biomass of each part of the plant; D is diameter at breast height; H is tree height; and a and b are constants. Finally, the biomass of the tree layer was calculated according to the regression equation.

Measurement of vegetation organic carbon indicators.
The harvested samples were crushed three times using the quadratic method, sieved using a 100 mesh, sampled, and dried to a constant weight in an oven at 80˚C. The carbon content (g�kg -1 ) was determined using the K 2 Cr 2 O 7 -H 2 SO 4 external heating method [31,32]. Three replicates of tree, shrub, and herb samples were analyzed for each lithology. The weighted average carbon content (g�kg -1 ) was calculated as the product of the carbon content of each part (or organ) and its weighting factor, that is, the ratio of the biomass of each part (or organ) to the overall biomass [10]. Carbon storage (g) was calculated as the product of the carbon content (g�kg -1 ) and biomass (kg), and convert carbon storage (t) to tons in every sample plot. Carbon density (t�hm -2 ) was calculated as the carbon storage divided by the area of the sample plot.

Determination of organic carbon in the litterfall accumulation.
A square sampling frame measuring 1 m 2 was prepared using a thin iron strip (1 mm thick and 10 cm wide), and 30 small randomly selected quadrats in each plot were sampled for the existing stock of litterfall. The square sampling frame was placed horizontally during sampling, and the existing stock of litterfall within the area of the frame, projected vertically in the sample plot, was collected as described elsewhere [18]. The litterfall was collected, brought to the laboratory, dried to a constant weight in an oven at 80˚C, and weighed. The litterfall accumulation was calculated based on dry weight. In addition, some dried samples were used to measure carbon content. Carbon density was calculated as the product of the litter fall accumulation per unit area (t�hm -2 ) of the quadrat and the mean carbon content of the litterfall in the corresponding quadrat [25,33].

Determination of soil organic carbon content.
In the present study, soil organic carbon content, organic carbon density, and organic carbon storage were measured using the soil volume actual measurement method in two layers (0-20 cm and >20 cm), as described elsewhere [12]. Following collection, the samples were placed into sealed plastic bags, brought to the laboratory, air dried, and subjected to soil organic carbon measurement. To determine bulk weight, samples were collected with a ring knife sampler [15]. Organic carbon content was determined using the K 2 Cr 2 O 7 -H 2 SO 4 external heating method [31]. Gravel content was determined by passing the sample through a 2 mm soil sieve [15]. Soil organic carbon storage (400 m 2 sample plots) was calculated using the following equation: where i represents the two soil layers (0-20 and >20 cm); CS i is the organic carbon storage of the corresponding soil layer; C i is the soil organic carbon content of the corresponding soil layer (%); D i is the bulk weight of the corresponding soil layer (g�m -3 ); G i is the gravel content (>2 mm) of the corresponding soil layer (%); and V i is the measured soil volume in the corresponding soil layer (m 3 ). Carbon density was calculated as the carbon storage divided by the sample area.

Calculation of indicators related to organic carbon in ecosystems.
Organic carbon storage in the ecosystem was calculated as the sum of carbon storage in vegetation, litterfall, and soil in the sample plots. Organic carbon density in the ecosystem was calculated as the ratio of organic carbon storage to sample plot area. Net organic carbon density was calculated as the difference between organic carbon density in the ecosystem and bare ground. The carbon sequestration rate of vegetation was calculated as the organic carbon density of vegetation divided by the recovery time. The carbon sequestration rate of the ecosystem was calculated as the difference between the organic carbon density of ecosystem and that of bare ground divided by the recovery time.

Data processing
The data were statistically analyzed using Excel and SPSS 17.0. One-way ANOVA and the least significant difference (LSD) test were used for comparisons of different datasets. The significance level was set at α = 0.05.

Vegetation biomass of the two plantation under different lithologies in karst regions
3.1.1 Regression models for relative growth relationships of vegetation tree layers under different lithologies. As shown in Table 2, the relationship of the biomass of each part of the tree (W) with DBH (D) and tree height (H) could fit well to the equation W = a

PLOS ONE
Correlations with this equation reached the level of significance at p<0.05. Dolomitic sandstone showed the optimal goodness-of-fit, followed by limestone, while dolomite showed the worst goodness-of-fit. These results indicate that plant growth was more consistent on dolomitic sandstone and more irregular on dolomite. Hence, dolomitic sandstone habitat is relatively consistent, followed by limestone habitat, while the dolomite habitat is complex and variable. Table 3, the overall vegetation biomass of B. cinnamon of the same or similar age was high in dolomite habitat, moderate in limestone habitat, and low in dolomitic sandstone habitat. Meanwhile, the overall vegetation biomass of C. funebris of the same age was higher in limestone than in dolomite habitat, indicating differences in species adaptability to different lithologies. The ratio of biomass of each layer to the vegetation biomass was high in the tree layer and low in the grass and shrub layer. The percentage of tree layer biomass to vegetation biomass was as follows: B. cinnamon in dolomite = 85.01%, C. funebris in dolomite = 76.37%, B. cinnamon in dolomitic sandstone = 50.08%, C. funebris in dolomitic sandstone = 88.25%, B. cinnamon in limestone = 71.10%, and C. funebris in limestone = 79.51%. All values were above 50%, indicating that vegetation biomass was dominated by the tree layer. At the same age, the biomass of B. cinnamon forest was greater than that of C. funebris forest in dolomite habitat, while the opposite result was noted in limestone habitat, indicating that B. cinnamon is highly adapted to dolomite, whereas C. funebris is highly adapted to limestone. Table 4, the overall ratio of aboveground biomass to belowground biomass was in the order of tree layer>shrub layer>grass layer, indicating that grass and shrubs show a greater adaptability to karst habitats. The aboveground biomass to belowground biomass ratio of trees and vegetation was high in dolomite, moderate in dolomitic sandstone, and low in limestone. Therefore, dolomite has better quality of habitats. The ratio of aboveground biomass to belowground biomass of C. funebris in the tree layer on limestone was the highest at 3.78. The plantation plot was a demonstration site of C. funebris blast plantation in a rocky desertified habitat. In this type of plantation restoration, a planting hole is created via micro-blasting on the rock, and the planting soil is collected from thin soil on the bare ground. The planting hole is small, and the surrounding rocks restricts root growth. Therefore, the belowground rocky habitat has a greater influence on root growth, and the belowground habitat itself is influenced by lithology. These results indicate the selectivity of species to lithology, which is a challenge during plantation restoration.

Carbon content of vegetation in 2 plantation forests under different lithologies.
As shown in Table 5, the carbon content of vegetation was in the order of tree layer>shrub layer>herb layer, and the carbon content of each tree part was in the order of stem>bran-ches>leaves>roots (tree layer). There was no significant difference between the aboveground and belowground carbon content of the shrub layer, indicating that lithology did not have a great impact on carbon content in forests. The weighted average carbon content of the stem and root of C. funebris in the tree layer at 38 years of age was significantly higher than that at 16 years of age, indicating that age had a great impact on carbon content. Table 6, the carbon density of vegetation in B. cinnamon forests of the same age was in the order of dolomite>limestone>dolomitic sandstone. Meanwhile, in C. funebris forests of the same age, the carbon density of vegetation was higher in limestone than in dolomite, indicating that lithology had a greater influence on the vegetation carbon pool in this forest. In C. funebris forests of different ages, carbon density was higher in older forests, indicating that age has a greater impact on vegetation carbon sink. The carbon density was higher in the tree layer, followed by the shrub and herb layers. Table 7, the biomass of litterfall in B. cinnamon forest was in the order of dolomite>limestone>dolomitic sandstone. Meanwhile, in C. funebris forest, the biomass of litterfall was higher in limestone and dolomitic sandstone than in dolomite. The carbon content of litterfall was higher in limestone and lower in dolomite and dolomitic sandstone. The carbon density was the highest in dolomitic B. cinnamon forest and the lowest in dolomitic C. funebris forest. Therefore, different lithologies significantly affected the litterfall accumulation and its carbon density.

Characteristics of the soil organic carbon pool in the bare ground under different lithologies in karst region
3.4.1 Soil volume weight, gravel content, volume, and organic carbon content of bare ground under different lithologies in karst region. As shown in Table 8, the soil volume and gravel content of bare ground tended to increase with the deepening of soil layer. In contrast, the soil volume and soil organic carbon content displayed a decreasing trend. The soil volume weight of bare ground was high in dolomite, moderate in dolomitic sandstone, and low in limestone. The gravel content of bare ground was high in dolomitic sandstone, moderate in limestone, and low in dolomite. The soil volume of bare ground was high in limestone and dolomitic sandstone and low in dolomite. The soil organic carbon content of bare ground was high in dolomite, moderate in limestone, and low in dolomite. Soil organic carbon content in the 0-20 cm layer was high in dolomite, followed by limestone, and low in dolomitic sandstone. Therefore, soil quality varies with lithology. Table 9, the carbon storage and density in the 0-20 cm soil layer of the bare ground were high in dolomite, moderate in limestone, and low in dolomitic sandstone. The carbon storage and density in the >20 cm soil layer were high in limestone, moderate in dolomite, and low in dolomitic sandstone. The carbon storage and density in the two layers together were consistent with the variation in the 0-20 cm soil layer. The carbon storage and density in the 0-20 cm soil layer were higher than those in the >20 cm soil layer. Therefore, soil carbon exhibits a surface clustering feature. Overall, carbon storage and density of dolomitic sandstone were lower, indicating the poor carbon sequestration capacity of soil developed from this parent material. Table 10, the soil organic carbon content and volume were higher in the 0-20 cm soil layer than that in the >20 cm soil layer under the three lithologies of the two plantations, while soil volume weight and gravel content followed the opposite trends. Organic carbon content was generally higher in the C. funebris forest than in the B. cinnamon forest under the same lithology. Organic carbon content of the 0-20 cm soil layer was higher in dolomite and limestone and lower in dolomitic sandstone. Organic matter in this layer is favorable for seedlings at the early stages of artificial vegetation restoration in dolomitic and limestone habitats. The overall soil volume was slightly higher in the B. cinnamon forest than in the C. funebris forest under the same lithology, while gravel content followed the opposite trend. Table 11, the storage and density soil organic carbon of each stand were higher in limestone, followed by dolomite, and lower in dolomitic sandstone. The storage and density of soil organic carbon in the C. funebris forest were higher than those in the B. cinnamon forest on dolomitic sandstone and limestone. No evident difference in the storage and density of soil organic carbon between the two vegetations was observed on dolomite. Therefore, dolomitic sandstone and limestone are more favorable for C. funebris growth. Under the same lithology, the storage and density of soil organic carbon in the 0-20 cm soil layer were significantly higher than those in the >20 cm soil layer, further indicating that soil organic carbon shows a surface clustering feature.

Carbon density and carbon sequestration rate of the two plantation forest ecosystems under different lithologies in karst regions
As shown in Table 12, under the three lithologies, among the components of ecosystem, carbon density was high in vegetation, followed by soil, and low in litterfall. Vegetation carbon

Gravel content (%) Soil volume (m 3 ) Soil organic carbon content (g�kg -1 ) (g�kg -1 )
0-20 cm >20 cm 0-20 cm >20 cm 0-20 cm >20 cm 0-20 cm >20 cm sink occupied the dominate position, indicating that vegetation is the most effective or contributes the most to the carbon sink in karst forest restoration. The net carbon density of the ecosystem in B. cinnamon forest of the same or similar age was high in dolomite, followed by limestone, and low in dolomitic sandstone. Meanwhile, in C. funebris of the same age, the net carbon density of the ecosystem was higher in limestone and lower in dolomitic sandstone. Therefore, dolomitic sandstone plantation forests have a poor carbon sequestration capacity. The carbon sequestration rate of the vegetation and ecosystem was high in dolomite, moderate in limestone, and low in dolomitic sandstone in B. cinnamon forest, whereas high in limestone, moderate in dolomite, and low in dolomitic sandstone in C. funebris forest. These trends support the poor carbon sequestration capacity of dolomitic sandstone plantation forests. Furthermore, the carbon sequestration capacity of plantation forest ecosystems differed with lithology, species, and plantation age, reflecting the complexity and diversity of karst lithological habitats and the difficulty of ecological restoration and carbon sink assessment in these regions.

Vegetation biomass characteristics of the two plantation forest ecosystems under different lithologies in karst region
The biomass of herbs, shrubs, and standard trees was determined using the harvesting method. However, the belowground biomass of trees was difficult to obtain because of the thin soil layer in the karst region, continuous distribution of rock in the belowground habitat, and general distribution of roots in rock pores and fractures. To overcome the difficulties of harvesting and to achieve accurate biomass determination, the present study used regression simulations of standard trees biomass in the harvested sample plots, as described in previous studies in the same area [10,[27][28][29][30]. We used the developed regression equation to determine the biomass of the tree layers. To verify the accuracy of the regression equation, the biomass of the tree layer was measured using the average standard tree method [24] and compared with the value obtained using the regression method. The results were consistent between the two methods, indicating the reliability of results obtained in the present study. Therefore, to measure the biomass of B. cinnamon and C. funebris in the same or similar study areas, the regression equations obtained in the present study can be used for accurate biomass determination based on surface dendrochronological indicators (i.e., DBH and H); this information can provide insights into vegetation biomass, particularly under different lithologies. Native forests in karst regions have a low biomass [10,27]. The biomass of plantation forests estimated in the present study was generally consistent with the reported values. For instance, the biomass of tree layer in B. cinnamon forests at 13-16 years of age ranged from 6.00 to 24.11 t�hm -2 , which is consistent with the values for artificial camphor forests at 16 years of age (tree layer: 15.18 t�hm -2 , vegetation: 17.21 t�hm -2 ) in Guiyang in the same karst region [9]. The biomass of the 16-year-old B. cinnamon plantation forest in the present study was also similar to that of the native karst shrubland (34.28 t�hm -2 ) [10,34] over the same restoration time but was lower than that of the 16-year-old plantation forest of camphor (9.24-24.91 t�hm -2 ) within the same plant family and in the same latitude-normative landform [35]. Furthermore, the biomass of 13-16-year-old B. cinnamon vegetation in the present study (11.98-28.36 t�hm -2 ) was much lower than that of the 18-year-old plantation forest of camphor within the same family and in the same latitude-normative landform (111.08 t�hm -2 ) [36]. These results indicate that plantation forests in karst regions produce lower than naturally restored native forests and plantation forests in a normal landform.
As shown in Table 13, the ratio of biomass of each tree part to the total biomass of the tree layer was in the order of stem (31.51%-41.63%) > branch (18.53%-30.06%) > root (18.52%-28.32%) > leaf (10.11%-20.59%). Therefore, stem and branches contributed the most to the overall biomass, although the contribution of roots cannot be ignored. The proportion of leaves was higher in C. funebris (16.40%-20.59%) than in B. cinnamon (10.11%-13.92%), while the branches followed the opposite trend. The proportion of leaves in B. cinnamon was higher

PLOS ONE
than that in 18-year-old camphor within the same family and genus and in a normal landform at the same latitude (4.70%) [36]. Furthermore, these values were higher than the proportion of leaves in 7-year-old camphor (7.24%-10.78%) in a normal landform at the same latitude [35] but similar to that in 16-year-old camphor (13.39%) in Guiyang within the same karst region [9]. Moreover, the weight ratio of B. cinnamon and C. funebris leaves was consistent with the ratio of leaf biomass to total biomass in the same karst region (12.28% for rowan, 8.24% for pepper, 43.04% for cypress, and 23.89% for Dodonaea viscosa) [5], suggesting higher leaf weight ratio in karst regions. However, the generalizability of this conclusion must be tested using various methods in future research. The ratio of aboveground biomass to belowground biomass reflects the growth strategy of plants, and it is to one of the core aspects of plant life history response theory [37]. The ratio of aboveground biomass (including leaves and stems) to belowground biomass of native plants is higher in karst region than that in normal landscapes [10], and our results are consistent with this conclusion. Here, the ratio of aboveground biomass (including leaves and stems) to belowground biomass ranged from 1.42 to 1.89 for the herb layer, 1.53 to 2.29 for the shrub layer, 2.97 to 3.78 for the tree layer, and 2.52 to 3.15 for the total vegetation. However, these values are lower than the global forest stem biomass-to-root biomass ratio of 3.13 to 5.88 and the Chinese forest stem biomass-to-root biomass ratio of 2.78 to 5.88 [38], indicating that the root biomass of karst plantation forest ecosystems accounts for a larger proportion of the overall biomass [10]. In addition, to adapt to the harsh ecological environment, plants in karst regions adopt an adaptive strategy of expanding root biomass to obtain more nutrients and groundwater.

Carbon sequestration characteristics of the two plantation forest ecosystems under different lithologies in karst regions and insights
In the present study, the carbon density in the litterfall accumulation was 0.44 to 1.29 t�hm -2 , being lower than the values from the herb stages to climax stages layers (1.53 to 4.97 t�hm -2 ) in naturally restored karst forests in the same region [18]. Therefore, carbon storage in the litterfall accumulation of an artificially restored karst forest was much lower than that of a naturally restored forest. The organic carbon content of the litterfall accumulation was 38.33% to 43.81%, which was comparable to or slightly higher than that of the native forest litterfall (37.78%-39.93%). Therefore, litterfall biomass in plantation forests (1.15-3.27 t�hm -2 ) is lower than that in native forests (3.83-13.16 t�hm -2 ) due to the simple species sources of litterfall in the former. In plantation forests, understory grasses and shrubs are artificially removed during the nurturing process, while in native forests, litterfall has abundant species and stable understory grasses and shrubs. Therefore, in karst regions, due to the harshness of the surface ecology, should the grass-shrub system in the understory of the dominant species be removed entirely or retained partially at the nurturing stage to create good hydrothermal conditions in the understory when plantation forests are cultivated? How much understory should be retained? It is worthwhile to study these issues in depth in future.
To objectively reflect the characteristics of the soil organic carbon pool [10], soil organic carbon content, soil volume weight, gravel content, and volume were actually measured. The organic carbon content of the 0-20 cm soil layer in the studied karst forest (32.37-39.88 g�kg -1 ) was similar to that reported in the Maolan Nature Reserve (20.89-203.5 g�kg -1 ) range [12] but higher than that reported in the Huajiang River (25˚38 0 N, 105˚38 0 E) [39] (12.68 g�kg -1 in grassland, 18.06 g�kg -1 in shrubland, and 20.96 g�kg -1 in forest) as well as the minimum content reported in the karst region of Southwest China (7.3 g�kg -1 ) [40]. Therefore, soil organic carbon content in karst regions is more variable. Hence, to assess soil organic carbon storage in karst regions, the actual measurement method is appropriate; otherwise, the measured results may be inaccurate. Therefore, the actual measurement has important methodological significance in the assessment of ecosystem carbon in karst regions. Further, the density of soil organic carbon in the present study ranged from 1.55 to 3.01 t�hm -2 , being lower than that reported in naturally restored forests of Mao Lan (8.26-18.80 t�hm -2 ) [12] and much lower than that of the camphor plantations of the same age in Guiyang (115.37 t�hm -2 ) [9]. Therefore, the soil organic carbon content in restored plantation forests in areas with intense rocky desertification is low. Simultaneously, compared with that in the bare ground, the soil carbon density in forests only increased by 0.16 to 1.28 t�hm -2 , indicating no significant increase under each lithology. Therefore, the plantation soil has a poor carbon sequestration capacity due to low soil volume and high gravel content in areas with intense rocky desertification in karst regions. Furthermore, the lower volume of soil must provide nutrients to plants, leading to higher soil respiration and organic carbon activity in forests of karst regions than in those in forests of normal landscapes [12]. Hence, soil carbon sequestration is poor. These factors reflect the limitations and unique features of karst plantation forests to increase carbon sinks.
The aboveground carbon content of the herbaceous layer in the present study ranged from 37.71% to 39.56%, being higher than that in naturally restored forests (37.13%) [41]. The carbon content of the shrub layer in the present study ranged from 44.64% to 47.11%, being close to that in natural forests (45.93%) [41]. The carbon content of various parts of the tree layer in the present study ranged from 44.44% to 51.22%, being different from the value commonly used in large-scale studies (45% or 50%). Therefore, the carbon content of plants in karst areas is highly variable. Furthermore, carbon content is one of the important factors affecting plant carbon storage; thus, a fixed carbon content should not be used in the estimation of vegetation carbon storage in karst regions, and the actual measurement method should be adopted instead. In the present study, the vegetation carbon density in young B. cinnamon forests ranged from 5.47 to 13.31 t�hm -2 , being consistent with that in artificial camphor forests of same age in the same area (9.98 t�hm -2 ) [9]. The vegetation carbon density in young and middleaged C. funebris forests in the present study was 9.64-11.29 t�hm -2 and 15.89 t�hm -2 , respectively. In a previous study [42], the average carbon density in young and middle-aged C. funebris forests in Guanling, Zhenfeng, and Huajiang in the same karst area was respectively 17.34 and 22.00 t�hm -2 ; these value are higher than those obtained in the present study. This may be because before restoration, the vegetation sample plots in the present study were in an intensely rocky desertified habitat, where plant growth is difficult. However, the carbon density of vegetation in plantation forests in the present study and that reported in the study by Ma (2020) [42] were lower than the values in naturally restored shrublands (25.60 t�hm -2 ) of the same age [10,34]. Therefore, habitat, tree species, forest age, and restoration methods in karst regions affect vegetation carbon density. In addition, the carbon density of forest vegetation in karst regions shows extremely spatial variability, leading to complexity and uncertainty in the assessment of carbon pools in these ecosystems. Ecosystem carbon density constitutes the carbon density of litterfall, soil, and vegetation. The ratio of vegetation carbon density to ecosystem carbon density was 79.94% and 79.65% for dolomitic B. cinnamon and C. funebris forests, 70.49% and 87.40% for dolomitic sandstone B. cinnamon and C. funebris forests, and 74.13% and 75.67% for limestone B. cinnamon and C. funebris forests, respectively. Therefore, the ecosystem carbon pool of artificially restored plantations in karst areas with intense rocky desertification is dominated by vegetation, as opposed to the dominance of soil in the ecosystem carbon pool in restored plantations on normal landscapes. This is mainly because normal landscapes have a large soil volume with rich organic carbon pools. Conversely, soil volume in the intensely rocky desertified areas is small, and the weathering of rocks into soil in such areas is a long process. Hence, it is difficult to achieve a significant increase in soil volume within the vegetation restoration period. The ratio of vegetation carbon storage to soil carbon storage in the present study ranged from 3.53 to 9.03, being higher than the average global value of 0.46 [43] and Chinese average of 0.36 [33]. This further confirms the greatest contribution of vegetation and the smallest of soil to the carbon pool of forest ecosystems in karst regions. These results have important implications for the establishment of carbon sink forests in karst regions. As such, it is necessary to focus on vegetation carbon sinks and pay particular attention to the selection of species and their cultivation techniques, in addition to further research and development.
The carbon density of ecosystem in the present study ranged from 7.76 to 16.65 t�hm -2 , being lower than that of evergreen and deciduous broad-leaved forests on a normal landscape in the same zone [44][45][46]. Therefore, plantation forests in rocky desertified areas are low-carbon ecosystems, with carbon density lower than the global average for forest vegetation (86.00 t�hm -2 ). In the present, the net carbon density of different lithological ecosystems was high in dolomite (14.83 t�hm -2 ), followed by limestone (9.47 t�hm -2 ) and dolomitic sandstone (6.36 t�hm -2 ). In B. cinnamon forests of the same age, lithology had a great influence on ecosystem carbon sink, possibly due to the differences in the adaptability of plants under various belowground lithological conditions [2]. This finding has important implications in the establishment of carbon sink forests in karst areas. Specifically, the strategy of "suitable lithologies and trees" should be adopted when selecting tree species for carbon sink forests. The carbon sequestration rate of the two plantation ecosystems in the present study (0.44-0.93 t�hm -2 �a -1 ) was rather different from the previously reported values in the same region (0.08 t�hm -2 �a -1 for cypress, 0.69 t�hm -2 �a -1 for Dodonaea viscosa, 0.28 t�hm -2 �a -1 for pepper, and 1.18 t�hm -2 �a -1 for rowan) [5], indicating that lithology, species, and habitat in karst areas significantly affect ecosystem carbon sequestration. The vegetation carbon sequestration rate in the present study (0.42-0.83 t�hm -2 �a -1 ) was slightly lower than that at all stages of natural recovery (0.45-1.70 t�hm -2 �a -1 ) [10] but significantly lower than that of 5-14-year-old alder plantations in karst areas and natural subtropical forests in Huayuan County, Xiangxi (3.68-5.33 t�hm -2 �a -1 ) [47]. Moreover, the carbon sequestration rate of tree layer in the present study (0.22-0.72 t�hm -2 �a -1 ) was lower than that in a 7-year-old camphor forest (0.88-1.60 t�hm -2 �a -1 ) in Jiangxi [35] (285 0 0 46@N, 117˚33 0 32@E) with the same silvicultural density (1000-2000 plants�hm -2 ). Therefore, the carbon sequestration capacity of karst ecosystems is lower than that of normal ecosystems due to the low biomass of vegetation and the low volume of soil in the former. In summary, although the carbon sequestration capacity of karst plantation forest soil is much lower than that of vegetation, lithology and its soil influence plant growth. Specifically, as karst regions have a small volume of soil, plant roots mostly use rocks as the growth space. Furthermore, the internal spatial morphological structure of rocks and its nutrient and water endowment patterns are different, all of which strongly affect plant physiological and ecological functions; thus, different lithologies lead to different plant growth characteristics [2], affecting the carbon sequestration capacity of vegetation. Nonetheless, the carbon density in the top community ecosystem in karst areas (81.31 t�hm -2 ) was higher than the average carbon density of evergreen and evergreen deciduous broad-leaved forests in China (73.68 t�hm -2 ) but close to the global forest vegetation carbon density (86.00 t�hm -2 ) [10]. Hence, karst forests can have a high carbon sequestration potential as long as they can reach the top community, and their large area produces a significant carbon sink effect of vegetation. However, the key question is how to build a carbon sink forest? The present study demonstrates that lithology is the major factor affecting vegetation carbon sink and belowground rocks and their soils have a great influence on plant growth and development. Therefore, the strategy of "suitable lithologies and trees" is of great importance to enhance carbon pools in karst areas to establish carbon sink forests. In addition, vegetation cultivation techniques tailored to different lithology should be developed.

Conclusion
Actual measurements and statistics based on soil volume, soil gravel content, plant carbon content, and biomass in karst areas are fundamental to ensure the accuracy of carbon sink assessments in forest ecosystems. In particular, fixed values of plant carbon content should not be used in carbon sink assessments. Karst plantation forests are the ecosystem with low biomass, carbon density, and carbon sequestration. To cope with the harsh karst environments and obtain more nutrients and water, plants have adopted a strategy to expand the proportion of their root and leaf biomass; thus, the ratios of plant root and leaf biomass to the total biomass are higher in karst regions than in normal landscapes. At the nurturing stage, the strategy of the removal or retention of the understory grass-shrub system in karst plantation forests should be different from that in plantation forests in normal landscapes, and these strategies warrant further research. In restored plantation forest ecosystems in intensely rocky desertified karst areas, vegetation is the major contributor to carbon sinks, and the contribution of soil is small. The carbon sequestration capacity of karst plantation forest ecosystems is not as good as that of native forests and forest ecosystems in normal landform; nonetheless, the carbon sink effect of these karst ecosystems remains significant if they eventually evolve toward being the top community. Lithology, forest age, and forest structure have a great influence on the carbon pool characteristics of karst plantation forest ecosystems. Therefore, the strategy of "suitable lithologies and trees" should be adopted during the establishment of carbon sink forests in karst regions.